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Abstract The fully differential angular distribution for 
the rare flavor- changing neutral current decay B° — > 
K*° (— > K~ix + ) /i + is studied. The emphasis is placed 
on accurate treatment of the contribution from the pro- 
cesses B® — > K*°(— > K~ 7r + ) V with intermediate vec- 
tor resonances V = p(770), w(782), </>(1020), J/^, ^(25), 
. . . decaying into the pair. The dilepton invariant- 

mass dependence of the branching ratio, longitudinal 
polarization fraction /j, of the K*° meson, and forward- 
backward asymmetry Afb is calculated and compared 
with data from Belle, CDF and LHCb. It is shown that 
inclusion of the resonance contribution may consider- 
ably modify the branching ratio, calculated in the SM 
without resonances, even in the invariant-mass region 
far from the so-called charmonia cuts applied in the ex- 
perimental analyses. This conclusion crucially depends 
on values of the unknown phases of the B° —> K*° J/ip 
and B° —> K*° ip(2S) decay amplitudes with zero he- 
licity. 

PACS 13.20.He; • 13.25.Hw; • 12.40.Vv 



1 Introduction 

The investigation of the rare decay 
B { j^K*° (-> K- 7T+) fi+ n~ 



induced by the flavor-changing neutral current (FCNC) 
transition b — > s/i + /i _ is an important test of the stan- 
dard model (SM) and its extensions (see [1] for a re- 
view) . The phenomenology of this decay mode has been 
discussed by many authors, e.g. see Refs. [i-27|. 



This decay takes place in a very wide region of 
dimuon invariant mass squared, q 2 = (q++q-) 2 , namely 
4rn 2 < q 2 < <j 2 lax = {ms — m K ') 2 ■ The light vec- 
tor resonances p(770), w(782), </>(1020) (and their ra- 
dial excitations) and the cc resonances J/tp, ij)(2S) (and 
higher states) are also located in this region. Thus at 
q 2 fa rriy the decay B® — > K*° /i + /z~ can go through 
the hadronic weak decay B® — > K*° V, followed by 
the dimuonic annihilation of vector meson V. All reso- 
nances with my < tub — TriK* make a contribution to 
this mechanism. Therefore both the nonresonant and 
resonant parts can contribute to the total amplitude of 
the decay B° -> K*° ^+ fi~. 

Main attention in literature has been paid to de- 
scription of the nonresonant amplitude of the decay 
B° d -> K*° m + (J>~~ in the region 1 GeV 2 < q 2 < 6 GeV 2 . 
In this region, using the QCD factorization (QCDf ) 28, 
[29| , one can perform a systematical calculation of non- 
factorizable corrections to "naive factorization approxi- 
mation" (NFA) and spectator effects 3(3, 3l[ . At larger 
dimuon masses, at about q 2 > 14 GeV 2 , the QCDf 
and the light-cone sum rules (LCSR) methods are not 
applicable. For the estimation of non-factorizable cor- 
rections, an operator product expansion in powers of 

< 



can be used [32], [33j. In the region 10 GeV 



g 2 < 13 GeV 2 the non-factorizable effects due to 
soft-gluon emission have been included in [34] . 

Often the resonant contribution to amplitudes of 
rare decays of B-meson is modeled in terms of the Breit- 
Wigner functions for the resonances 35-3^|. In these 
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references the resonance corrections are added to the 
perturbative loops of charm quarks. Note an original 
approach of Ref. 40] for the inclusive B — > X s 
process, in which dispersion relation exploiting exper- 
imentally measured cross section a(e + e~ — > hadrons) 
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has been applied to account for the resonance terms 
(see also |41|). 

These approaches more often apply to the inclu- 
sive decays B — > X s £ + £~ and use information on the 
B — s- X s J/tp and B —> X s tp(2S) branching ratios for 
description of the resonant contribution. Sometimes, 
such approaches are extended to the exclusive decays 
B K{K*)i+£- @, [H HI, in which the branch- 
ing ratios for exclusive decays B — > K(K *) J/tp and 
B — > K(K*) tp(2S) are used. In these studies, carried 
out in framework of the NFA, additional factors ky 
are introduced into the resonant terms to adjust the 
branching ratios for the decays, for instance, 

BR(B -S- K*V -> K* £+£-) 

= BR (B ->■ K* V) BR (V ->■ l+T) , 

where the right-hand side is taken from experiment. 

Recall that, in general, the process B — > K* V is 
characterized not only by the branching ratio. The de- 
cay of a Bj meson into a pair of vector mesons, B® 
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K*° V, is described by three complex amplitudes 
In the transverse basis 4^, 44 j , these decay amplitudes 
correspond to linearly polarized states of vector mesons, 
which are polarized either longitudinally (0) or trans- 
versely to the direction of their motion, being polar- 
ized in parallel (||) or perpendicular (_L) to each an- 
other. Overall, six real parameters describe three com- 
plex amplitudes Aq , Aj , and A\ . They could be chosen 
to be, for example, the branching ratio, l-Ag^ 2 , l^l 2 , 
axg(A^ /AX), arg(^Y l A o) and one overall phase arg(A^ 
The phase convention is arbitrary for an isolated decay 
B^ — > K*° V. Sometimes, this phase is chosen zero, 
arg(Aj') = 0. However, for certain B decays, this phase 
can produce meaningful and observable effects, such as 
for B — !• V Kj with J = 0,1,2,... . 

For example, in the decay B^ — > K*° <fi, the phase 
of the amplitude Aq has been measured with respect to 
the phase of the amplitude Aq of the decay B® —> 
A O *(143O)°0 and is equal to axg(A$ /A$) = 2.82 ± 



0.15 ± 0.09 [45l. l46j]. For the other vector resonances V, 



the corresponding relative phase has not been measured 
so far. 

At present in decay modes to the light resonances, 



higher states, ^(3770), </>(4040), -0(4160), ^(4415), the 
experimental information is absent. 

In the present paper for description of resonant con- 
tribution to the four-body decay B® — > K~ ir + fi + 
the available information on the helicity amplitudes for 
B — » K* V is used. The fully differential angular dis- 
tribution over the three angles and dimuon invariant 
mass is analyzed in the whole region 4m 2 < q 2 < 
(tub — nix*) 2 - The amplitude of this decay consists of 
the nonresonant amplitude in the SM model and the 
resonant amplitude. For the first amplitude we use the 
NFA, in which hadronic matrix elements are parame- 
terized in terms of form factors and the Wilson 
coefficients are taken in the next-to-next-to-leading or- 
der (NNLO) approximation. 

The resonant amplitude is expressed in terms of the 
invariant amplitudes SY% 3 for the decays B® — > K*° V. 
The information on the latter is taken from experiment 
if available, or from theoretical estimations. As men- 
tioned above, the phase 8$ = arg(^4^) for an isolated 
decay B^ — > K*° V is arbitrary. This phase may pro- 
duce observable effects in the decay B® —> K*° /j+ pT 
via the interference with the nonresonant amplitude. 
We investigate influence of the phases Sq for each reso- 
nance V on the differential branching ratio, longitudinal 
polarization fraction /l of K* and forward-backward 
asymmetry Apn- 

We also study two aspects of the resonant ampli- 
tude. The first one is related to the fact that the vector 
mesons are off their mass shells, therefore an off-mass- 
shell extension of the on-mass-shell B® — > K*° V am- 
plitudes is proposed. The second one is the choice of 
the vector-meson-dominance (VMD) model which de- 
scribes the transition vertex V7. We use two versions 
of the VMD model (called further VMD1 and VMD2) 
which result in rather different V j vertices. In partic- 
ular, in the VMD2 model, the V 7 vertex is suppressed 
compared to the VMD1 vertex in the region q 2 <C my. 

Results of the present calculations are compared 
with the recent data from Belle, CDF and LHCb exper- 
iments. Usually in these experiments the cc resonance 
contributions are removed by putting cuts on the in- 



K*° p(770) and B^ — > K*° cj(782) only the variant dimuon mass near the resonance mass q 



B° d 

branching ratio and longitudinal polarization fraction 
of the K* meson are measured, while the decays to 
radial excitations, w(1420), p(1450), w(1650), p(1700), 
• • • , have not been observed. At the same time, all am- 
plitudes of the decay B® — > K*° 0(1020) are known 
from experiment, while there is no information on B 
decays to radial excitations of 0(1020). For the B^ — > 
K*° J/tP and B° -> A"* iP(2S) decays, the full angular 
analysis has been performed. As for the decays to the 
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v- 

This assumption is used in the analyses of all ongoing 
and planned experiments. 

The paper is organized as follows. In Section 12.11 
the fully differential angular distribution is discussed. 
Nonresonant and resonant amplitudes in the transverse 
basis are specified in Section l2~2l Results for the depen- 
dence of observables on the invariant mass squared are 
presented in Section [3] Conclusions are drawn in Sec- 
tion [U In |Appendix~A"| calculation of the B° -t K*°V 
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Fig. 1 Definition of heiicity angies 9 M , 8k, and <j>, for the 
decay _B° ->■ K*° n~ . 



amplitudes for the off-mass-shell vector meson V is de- 
scribed. 



and the amplitude terms a k as 

ai =f3l\A \\a 2 =(3l(\A^ + \A^), 

-|A||| 2 ), a 4 = ^Im(A||Al), 



03 

a 5 = $ Re(A Ap , a 6 = /3 2 Im(iMI) , 

a 7 = Re(A llL A* ±L - A m A* ±R ) , 

a 8 = Rc(A 0L Al L - A 0R A* ±R ) , 

a 9 = lnv{Au L A\ L - A 0R A\ R ) , 

a 1Q = (l-f3l)(\A QL +A QR \ 2 + \A t \ 2 ), 

an = (1 - Pl)(\A llL + A m \ 2 + \A XL + A ±R \ 



where p^ = y 1 - 4m2/q 2 , m p 
mj is the mass of the B? meson, q 2 



is the mass of the muon, 
q 2 /m 2 



i 2 B , and 



2 Angular distributions and amplitudes for the 
K*° /Lt+ decay 



2.1 Differential decay rate 



The decay B° d Jf*V + M~, with K*° -> K-%+ on 
the mass shell 0, is completely described by four inde- 
pendent kinematic variables: the dimuon invariant mass 
squared, q 2 , and the three angles 8^, 9k, 4>- I n the he- 
iicity frame (Fig. [1]), the angle 6> M (6k) is defined as the 
angle between the directions of motion of /x + (if - ) in 
the 7* (K*°) rest frame and the 7* (K*°) in the B° 
rest frame. The azimuthal angle tj> is defined as the 
angle between the decay planes of 7* — > /i + /i~ and 
K*° — > K~ir + in the rest frame. The differential 
decay rate in these coordinates is given by 



AA* = A lL (q 2 )A* L (q 2 ) + A lR (q 2 )A* R (q 2 ) . 

Here i,j = (0, ||,_L), the a k dependent on products 
of the six transversity amplitudes A l( R }, and 
Aj_l(r), where L and R refer to the chirality of the lep- 
tonic current, as well as the seventh transversity ampli- 
tude A t . The latter amplitude is related to the time-like 
component of the virtual gauge boson, which does not 
contribute to the decay rate in the case of massless lep- 
tons and can be neglected if the lepton mass is small in 
comparison to the invariant-mass of the leptonic pair. 
Further, A = \(l,q 2 ,m 2 K ,) = (1 -q 2 ) 2 -2(l + q 2 )m 2 K , + 
rh K », rriK* = rriK* /iwb, where itik* is the mass of the 
K*° meson, and 



N=\V tb V*\ 



G F m 2 B a cm 



d 4 r 



dq 2 d cos 6 ad cos 0Kdcj) 



M m B N 2 q 2 VI— T a k {q 2 )g k (9^9 K A) , (1) 

04 7T * — ' 

k=l 



where the angular terms g k are defined as 

gi = 4 sin 2 M cos 2 6 K , 92 = (1 + cos 2 M ) sin 2 9 K 

93 = sin 2 9^ cos 2<f> sin 2 8k , 

94 = —2 sin 2 8^ sin 2 9k sin 2 <f> , 
.95 = — V2 sin 2 0^ sin 2 #x cos , 

,96 = — \/2 sin 2 # M sin 2 sin </> , 57 = 4cos( 
38 = — 2v / 2sin# /i sin2#K-cos0, 
(79 = — 2 V2 sin M sin 2 sin , 
g w = 2 cos 2 8 K ,gn= sin 2 6»if , 



,, sin 6»a' 



1 This means the narrow- width approximation for the K*° 
propagator: (k 2 — m 2 ^, + im K *r K »)~ x rj —in5(k 2 



'I- 



Here, 14,- are the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix elements 



49j , Gf is the Fermi coupling con- 
stant, a cm is the electromagnetic fine-structure con- 
stant. 

The longitudinal, parallel, and perpendicular partial 
widths are given, respectively, by 

dT L 



dq 



S2 



dr ± 
H4 1 



m B N 2 ^q 2 ^X[l3l\A \ 2 



q 2 



\A 



0L 



A 



0R\ 



3m 2 



= m B 



N 2 ^q 2 VX(f3l\A ± \ 



3m 2 



■ A 1 



(2) 



(3) 



(4) 
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The familiar muon-pair invariant-mass spectrum for — > 
K*° /i + /2~ decay can be recovered after integration over 
all angles as 

dr _ dr L dr l{ dr ± 

dq 2 dq 2 dq 2 dq 2 



(5) 



The fraction of K* meson polarization is [i — (L, || , _L)] 
» / 2 i d T{, .dr 

Mq)= w f w 

Integrating Eq. (JTJ over the variables cos # M and <j>, we 
obtain 

d 2 T dr /Z 

1^ Q~ = TF2 1 C0S " ^ 

dq d cos c'/f ag^ v 9 



+ J(Wl)(1-cos 2 ^)) 



(6) 



Integration of Eq. (JTJ over cos 6*#- and yields 



d 2 r 



dq 2 d cos p 



7V 2 /3 At g 2 yi|(2/3 2 |A | 2 si 



sin 2 0,. 



+^(| j4|| | 2 + |Ax| 2 ) (l W^) + ^(L4 t | 2 

+ |A0L + A B| 2 + |A|| L +A|| fl | 2 



-|Axi + A XR \ 



where dA^/dq 2 is the muon forward-backward asym- 
metry, 

d ^FB f , n s d 2 T J n 

— ~ ' sgn(cos^) ; —dcosu,, 



dA 



(p) 

FB 



dg 2 



(7) 



<7<j 



-,2 



M dq 2 d cos 9^ 



= m B N 2 plq 2 V 7 ! ^Rc (A,, L £ - A,, H ^ R ) , (8) 

and the normalized forward-backward asymmetry 
dA^/dq 2 is given as 



dA^l _ dA [ £l dr 



dq 2 dq 2 dq 2 



(9) 



Finally, the one-dimensional angular distribution in the 
angle <j> between the lepton and meson planes takes the 
form 

d 2 T 1 dr ( 1 , \ aI2 ) 

Wdj = ^w( 1 + 2 {1 - fL)A % cos2 * 

-Ai m sin 20), 



.4 



dr 



(2) _ 



dr ± 




j-(ii) 



dq 2 



m B 



N*fifVX\A ±m \*, 



(10) 

(11) 
(12) 



A 



,,„. = m B iV 2 /3f l( 7 2 yilm(A||Al)/^, 



(13) 



where the asymmetry A^ (q 2 ) is sensitive to new physics 
from right-handed currents beyond the standard model, 
and the amplitude Ai m (q 2 ) is sensitive to complex phases 
in the hadronic matrix elements. Sometimes A^\q 2 ) is 
called transverse asymmetry 



2.2 Resonant and nonresonant transverse amplitudes 

The effects of the long-distance contribution from the 
decays B° -> K*°V, where V = p°, w, (j), J/tp{lS), 
ip(2S), . . . mesons, followed by V — > fi + \i~ in the de- 
cay B% -> K*° fi + ^ are included through the VMD 
approach, as shown in Fig. [5J There is no unique way 



TC" 



Fig. 2 Nonresonant and resonant contributions to the decay 
amplitude. 



of introducing the V 7 transition, and one can use var- 
ious versions of VMD models which yield different V 7 
transition vertices. In one of VMD models (see [50j |. 
chapter 6) 



(7(9); V>\V(q); v) = -efvQvmv 9^ 



(14) 



where is the metric tensor, q is the photon (meson) 
four-momentum and Qy is the effective electric charge 
of the quarks in the meson V: 



1 

71' 



if) /jr J QtO 0/0' 



1 



3\/2 



'j/v — Qi/>(2S) — ■ ■ ■ — 



(15) 



The decay constant of the neutral vector meson fy can 
be extracted from electromagnetic decay width, using 



r 



A.-KO. 



V^re+e- 



3 ray 



em j?2 y^)2 
iyHy ■ 



(16) 



We will call this version VMD1. The vertex in Eq. (|TlJi 
follows from the transition Lagrangian 

C lV = -efyQymy A» . (17) 

Another model (called hereafter VMD2) originates from 

efvQv 



7 y 



2my 



. PC" y 



(18) 
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where V MV = <9 M K - d v V^ and = d»A u - d v A^ . An 
advantage of the Lagrangian (fT5|) is its explicit gauge- 
invariant form. Eq. (|18p gives rise to the V 7 transition 
vertex 



(7(g); n\V(q); u) 



efvQv 

TTiy 



(g V - q»qn , (19) 



The term oc q^q v does not contribute being contracted 
with the leptonic current eu{q-)^ ll v{q + ), and vertex 
(fT9|) is suppressed compared to (fl4|) at small q 2 <C my, 
i.e. in the region far from the vector-meson mass shell. 
Of course on the mass shell, q 2 — my, the VMD2 and 
VMD1 are equivalent. 

These two versions of VMD model have been dis- 



cussed earlier in Refs. [5lj, |52j. Note also that VMD2 
vertex follows from the Resonance Chiral Theory [53| 
and has been applied [54J when studying the reaction 
e + e~ — > 7r°7r°(77)7. 

Parameters of vector resonances are collected in Ta- 
ble EJ 

The nonresonant amplitudes are calculated in the 
NFA, with the short-distance NNLO Wilson coefficients, 
and nonperturbative B — ¥ K* transition form factors. 

Then the total amplitudes including nonresonant 
and resonant parts take the form 



.4 



1 



0L,R 



2 mp \J q £ \ v 
+2m h {Cf^ - C' 7 f) KQ (q 2 )) 
+8n 2 J2C V D v \q 2 )((l-q 2 



10.4 



i K , 



+A- 



(20) 



Ai 



L,R 



\C\\{q 2 )(ci v TC WA 



+ 2^(C^-C'f)^q 2 

v ) 



(21) 



A±l,r = V2A( CAq 2 )(c c 9V T c WA 



-2^(C% + Ctf)K ± tf: 
v / 



A, 



q 1 



(22) 
(23) 



where 

Dy(q 2 ) = q 2 - rhy + im v Py(q 2 ) 

is the usual Breit-Wigner function for the V meson res- 
onance shape with the energy-dependent width Iy(<? 2 ) 
[r v (q 2 ) = r v (q 2 )/m B ], my = m v /m B , f v = r v /m B , 
my(ry) is the mass (width) of a V meson and the form 
factors enter as 

C (q 2 ) = (l-q 2 - m 2 K ,)(l + niK^A^q 2 ) 
- A 2 (q 2 ) 



1 + rriK* 
C ll (q 2 ) = (l + m K ,)A 1 (q 2 ), 

1 + vriK* 

« (<Z 2 ) = ((1 -f + 3m 2 K ,)(l + m K ,)T 2 {q 2 ) 
-T 3 (q 2 ))((l-q 2 



(24) 
(25) 

(26) 



A 



1 - m K 

x(l + m K ,) 2 A 1 (q 2 )-XA 2 (q 2 ) 



m 2 K , 



K\\{q ) = 



K ± (q 2 



T 2 (q 2 



Ai(q 2 



-(1 - m K *), 



(27) 



(28) 



(29) 



In the above formulas the definition rhb = fni,(fi)/mB, 
m s = m s (n)/mB are used, and mj(^) [m s (/i)] is the 
running bottom (strange) quark mass in the MS scheme 
at the scale \i. 

The SM Wilson coefficients have been obtained in 
[lH at the scale /1 = 4.8 GeV to NNLO accuracy and 
equal 

CfiQi) = -0.304, C 9V (m)= 4.211, C 10A (//) = -4.103, 

Cg V — Cgy+Y(q 2 ), where Y(q 2 ) is quark-loop function 
given in Ref. [30( . Note that in the framework of the SM 
m b (fi)C' 7 f = m s (/i)C^. 

Further, A (q 2 ) 7 A x (q 2 ), A 2 (q 2 ), V(q% T^q 2 ), T 2 (q 2 ), 
Ts(q 2 ) are the B — > K* transition form factors. In the 
numerical estimations, we use the form factors from the 
LCSR calculation 0. 

In Eqs. (E01)-(E2), SY {i = 1, 2, 3) are the invariant 
amplitudes of the decay B^ — > K*° V. These ampli- 
tudes are calculated in |Appcndix A The coefficients 
Cy in the resonant contribution are 



Cy = Qvm Y fv (VMD1) , C V 



Qvfv 
my 



(VMD2) 
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Table 1 Mass, total width, leptonic decay width and coupling fv of vector mesons [55 
shown) . 



(experimental uncertainties are not 



mv-(MeV) i>(MeV) T v 



.(keV) fv(MeV) 





775.49 


149.1 


7.04 


221.2 




782.65 


8.49 


0.60 


194.7 


4> 


1019.455 


4.26 


1.27 


228.6 


J/i> 


3096.916 


0.0929 


5.55 


416.4 


V>(2S) 


3686.09 


0.304 


2.35 


295.6 


^(3770) 


3772.92 


27.3 


0.265 


100.4 


V>(4040) 


4039 


80 


0.86 


187.2 


j/>(4160) 


4153 


103 


0.83 


186.5 


i/)(4415) 


4421 


62 


0.58 


160.8 



Table 2 The numerical input used in our analysis. 



\VtbV t * s \ = 0.04026 
ii = mj = 4.8 GeV 
m c = 1.4 GeV 
m b (jt) = 4. 14 GeV 
= 0.079 GeV 



Gf 



1.16637 x 10- 5 GeV" 2 
Qom = 1/137.036 
m B = 5.27950 GeV 
t s = 1.525 ps 
m K , = 0.89594 GeV 



The energy-dependent widths of light vector reso- 
nances p, w and </> are chosen as in Ref. The up- 
dated branching ratios for resonances decays to differ- 
ent channels are taken from [55[ . For the cc resonances 
J/ip, t()(2S), ... we take the constant widths. 

In order to calculate the resonant contribution to 
the amplitude of the B^ — > K*° p + pT decay, one has 
to know the amplitudes of the decays B® —> K*° p, 



R0 



K* 



m 



K* 



m 



K*°J/4>, B° d 



K*° V»(2S). The information on the B Q d -> K*° 0, B° -> 
K*°J/ip, B°-> K*°ip(2S) decays can be taken from 
experiment [55[ . For the li ght resonances p and ui we use 
the theoretical prediction [56( for the decay amplitudes. 
At the same time, we are not aware of a similar predic- 
tion for the higher cc resonances, such as ^(3770) and 
so on, and therefore do not include these resonances. 

The parameters of the model are indicated in Ta- 
ble 1 



Results of the calculation for the 
K*° fi~ decay 



In Figs. [3] and [4] we present results for the dependence 
of several observables in the B® — > K*° p, + p~ decay on 
the dimuon invariant mass squared. The interval of q 2 
is taken from 4to^ to q^ ax — (tub — niK') 2 ~ 19.22 
GeV 2 . 

In order to adequately compare results of calcula- 
tions with experiments from Belle (KEKB) [57[, CDF 
(Tevatron) Hi H3| and LHCb [g^] we show the binned 



Table 3 Experimental bins in q 2 used in the data analyses 
[6p| (the first column), [57rl58|| (the second column) and in 
our calculation. 



q 2 (GeV 2 ) 


q 2 (GeV 2 ) 


0.10 - 2.00 


0.00 - 2.00 


2.00 - 4.30 


2.00 - 4.30 


4.30 - 8.68 


4.30 - 8.68 


10.09 - 12.86 


10.09 - 12.86 


14.18 - 16.00 


14.18 - 16.00 


16.00 - 19.00 


16.00 - ql ax 



predictions, i.e. the calculated values averaged over the 
experimental bins. These bins are indicated in Table [3] 

The solid line demonstrate the prediction of the SM 
without resonances for the g 2 -dependence of all observ- 
ables; the horizontal bars are the corresponding binned 
results. 

The dashed rectangles in Figs. [3] and |4] show possi- 
ble contributions of the intermediate vector resonances 
in each g 2 -bin. To obtain these rectangles the phases 
Sq = arg(A^) of the B — > K*V decay amplitudes with 
the zero helicity for V — p, u>, 4>, J/if>, ip(2S) have been 
varied from to 2ir. The upper (lower) side of the rect- 
angles corresponds to the maximal (minimal) value of 
the branching ratio in each bin. The dominant depen- 
dence on the phase 5$ comes from the cc resonances 
J/ip and ip(2S), while the phase dependence for the 
p, uj and <f> is of minor importance. This is explained 
by the much bigger contribution of the cc resonances 
compared to the light ones. 

Predictions of the VMD1 and VMD2 models (Sec- 
tion [212]) differ at relatively small q 2 . This is due to the 
fact that at small q 2 the high- lying resonances J/ ip and 
ip(2S) are far off their mass shells and in these kine- 
matical conditions the VMD1 and VMD2 models give 
different results. On the contrary, in the high q 2 region, 
close to the mass shells for the cc resonances, the results 
with VMD1 and VMD2 do not differ. 
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q 2 (GeV 2 ) q 2 (GeV 2 ) 



Fig. 3 Differential branching ratio as a function of q 2 . Solid lines correspond to calculation in the SM without resonances, 
and the horizontal bars (red online) are the corresponding binned (averaged over the experimental bins) results. In the first 
bin calculation is performed for 0.1 GeV 2 < q 2 < 2.0 GeV 2 . The dashed rectangles (blue online) show contribution of the 
intermediate resonances: left (right) panel - in the VMD1 (VMD2) model. The uncertainty of the resonance contribution is 
related to the unknown phase <5q of the resonance amplitudes (see the text). The upper side of the resonance rectangle (left side) 
is located at 3.54 and is not shown. The form factors are taken from [47| . The data Belle (KEKB) [57|. CDF (Tevatron) [58l.l59(| 
and LHCb 60] are shown by the circles, filled boxes and triangles respectively (the horizontal error bars are not shown - they 
are given in Table EJ. 
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Fig. 4 Longitudinal polarization fraction of K* meson (upper panel) and forward-backward asymmetry Ape (lower pan el) 
as functions of q 2 . Due to the choice of reference frame in Fig. ^ the forward-backward asymmetry Ape in Refs. [57h60|| is 
related to asymmetry in Eq. ((9j via A-p-g, = — dA^l/dq 2 . The other notation is the same as in Fig. [3j 
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In the experimental analyses certain cuts in q 2 are 
applied in order to suppress the charmonia contribu- 
tions (the so-called charmonia veto). Correspondingly 
in Figs. [3] and |4] we also use these cuts. As it is seen 
from Fig. [3j the charmonia resonances may contribute 
to the differential branching far beyond their pole po- 
sitions and beyond the charmonia veto, although their 
contribution crucially depends on the values of the zero- 
helicity phases Sq^ and Sq^ 2 ^. The observed in the 
VMD1 model sizable contribution of the cc-resonances 
at small values of q 2 occurs at certain values of these 
phases. Of course there exist values of 8^ and 5^ 2S ^ 
at which resonance contribution is small. 

Comparing these results with results of Refs. (3.1381 
|39( we note that in these papers the amplitudes of the 
decays B — >• K* J/ip and B -» K* ip(2S) are calculated 
in framework of the NFA which implies zero phases. 
In contrast, our calculation is free from this limitation 
since the amplitudes are taken from experiment and 
thus the phases Sq for all resonances V can be chosen 
arbitrary and independent. 

Note that long-distance effects due to the nonfactor- 
izable soft-gluon emission from the c-quarks have been 
included in Ref. [34|. Authors have shown that these 
effects lead to a modification, depending on the K* po- 
larization, of the Wilson coefficient Covin). The form 
factors and the Wilson coefficients in [34[ are different 
from the corresponding quantities in our calculation. 
In framework of this approach we estimated the differ- 
ential branching ratio in the region 1.0 GeV 2 < q 2 < 
9.0 GeV 2 , using the Wilson coefficients and form fac- 
tors of our work. Results of this calculation along with 
results of the present work are shown in Table |H 



Table 4 Differential branching ratio averaged over the two 
experimental bins. Results in our approach with resonances 
(VMD1 and VMD2) show variation of the branching ratio 
depending on the phases 5X. Values in the last row are ob- 
tained following approach 13411 for the central values of all 
parameters in Eq. (70) from |34| without taking into account 
uncertainties. 





dB/dq 2 (10- 


" 7 GeV" 2 ) 


q 2 (GeV 2 ) 


2.00 - 4.30 


4.30 - 8.68 


SM without resonances 


0.44 


0.63 


VMD1 


0.30 - 0.74 


0.43 - 0.94 


VMD2 


0.38 - 0.50 


0.48 - 0.82 


Ref. [34] 


0.44 


0.60 



As is seen, results of both approaches are compara- 
ble, at least in these intervals of q 2 and for a certain 
choice of phases Sq. Besides, our calculation in the SM 
without resonances in framework of the NFA turns out 



to be relatively close to the calculation following ap- 
proach of Ref. 0. 

As expected, the observables which are the ratios 
of the bilinear combinations of the amplitudes, such 
as /l or Afb, are less sensitive to the phases of the 
resonance contribution (see Fig. . Moreover, the latter 
is generally small independently of 5q, apart from the 
region of relatively small q 2 < 4 GeV 2 . 

It should be noted that for the inclusive decay B — > 
X s e + e~ , in Ref. [6l[ soft-gluon emission from the charm 
loop has been considered and long-distance corrections 
of 0(Aq CD /iti 2 ) to the differential branching and forward- 
backward asymmetry are shown to be small in the re- 
gion away from the cc resonances. Although this inclu- 
sive process is amenable to clean theoretical description 
for 1 GeV 2 < q 2 < 6 GeV 2 and q 2 > 14.4 GeV 2 [gj, 
the experimental information on differential branching 
ratio is still rather limited 



63, 641. The situation is not 



likely to improve in the nearest future. At the same time 
for the exclusive decays B — > K* (— » Kir) £ + l" a high 
statistics is expected at LHCb, and this will allow for 
measurement of various observables in these processes. 
Therefore the exclusive decays have quite a good po- 
tential for tests of the SM and search for effects beyond 
the SM. 

4 Conclusions 

The rare FCNC decay B° -> A'* (-> K~ tt+) /j+ fi~ 
has been studied in the whole region of muon-antimuon 
invariant masses. We performed calculations of the dif- 
ferential branching ratio, polarization fraction //, of the 
K* meson and forward-backward asymmetry ^4fb for 
the four-body decay B^ — > K~~ n + /i + fi~ . Main em- 
phasis in our study is placed on contribution of the 
intermediate vector resonances in the process B^ — > 
K*°(^ K~ir+)V with V = p(770), w(782), 0(1020), 
J/ip, ip(2S), decaying into the [i + [i~ pair. Various as- 
pects of theoretical treatment of this long-distance con- 
tribution have been investigated. 

One aspect is the choice of VMD model describing 
the Vj transition. We applied two different versions 
called VMD1 and VMD 2. In particular, VMD2 model 
explicitly obeys gauge invariance and the corresponding 
V 7 vertex is suppressed compared to the V 7 vertex in 
the VMD1 model at small values of the photon invari- 
ant mass, far from the vector-meson mass shell. This 
turns out to be important when calculating the contri- 
bution of the high- lying J/ip and ip{2S) resonances to 
the branching ratio at q 2 <C my. 

Another aspect is the structure of the vertex B% — >• 
K*° V with the off-shell vector meson. We used an off- 
mass-shell extension of the helicity amplitudes which 
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describe production of the on-shell vector mesons V. 
For the latter amplitudes all the experimentally avail- 
able information is used, and otherwise theoretical pre- 
dictions. 

The calculations are compared with recent data 
60] for the g 2 -dependence of the differential branching 
ratio, longitudinal polarization fraction of the K* me- 
son and forward-backward asymmetry. In these anal- 
yses the cuts around the pole positions of the char- 
monia resonances are applied to suppress these contri- 
butions (charmonia veto). Our calculation shows that 
the intermediate resonances may considerably modify 
the branching ratio, calculated in the SM without reso- 
nances in framework of the NFA, even in the q 2 -regioii 
located far from these cuts. Main contribution comes 
from the cc resonances J/ip and ip(2S). This conclu- 
sion, however, crucially depends on values of the un- 



known zero-helicity phases S and Sq for the de- 
cays B° -> K*° J/ip and P° -> K*° ip(2S) respectively. 
In view of these results one should keep in mind that 
the resonance contribution may imitate effects of new 
physics, the search for which is one of the goals of cur- 
rent and planned experiments on BaBar, CDF, Belle 
and LHCb. 

As for the observables Jl and Afb, these are less 
sensitive to the phases of the resonance contribution 
and moreover the latter is generally small independently 
of Sq (apart from the region of relatively small q 2 < 4 
GeV 2 in the VMD1 model). Apparently data for these 
observables are more appropriate for adequate compar- 
ison with the prediction of the SM without resonances. 



Then the phase structure of Si, S2 and S3 is: 

S k = J2\Sk j \e iv »e iS ^ (k = 1, 2, 3) , (A.2) 



where <pij, tp2j> and tpsj are the CP- violating weak 
phases and dij, 82 j, and 83 j are the CP-conserving 
strong phases. 

Using CPT invariance, we can represent the ma- 
trix element for the charge-conjugate decay B(p) — s- 
V 1 (q,e 1 )V 2 (k,e 2 ) as 



TJ GFm B /- S 2 

\ m B 



V2 



m 



~1 c 2 1 



(A.3) 



where Si, §2, and S3 can be derived from Si, S 2 , and S3 
by reversing the sign of the CP-violating phase. Note 
that if the B — > V\ V 2 decay is invariant under the CP 
symmetry, then Si = Si, S 2 = S 2 , and S3 = S3. On 
the other hand, if all CP-conserving phases of invariant 
amplitudes are equal to zero, then Si = S[*, S 2 = S 2 , 
and S 3 = Sa*. 

The helicity amplitudes in terms of three invariant 
amplitudes, Si, S 2 , and S3 are: 



P.\ 



o ^2 

Si g MV H 2~P»P v 



(A.4) 



Appendix A: Amplitudes of B — >■ K*V decays 

An important ingredient of the resonant contribution 
is amplitude of the decay of B meson into two vector 
mesons, B(p) —> V%(q, ex) + V 2 (k,e 2 ), with on-mass- 
shell meson V 2 (k 2 = m 2 ) and off-mass-shell meson 
^1 (l 2 7^ m i)- For the case of two on-mass-shell final 
mesons one can write the amplitude in the form |42| 



M 



GFm% / S 2 

7= I Vckm I ( Si g^u H 2" Vy.Vv 

V 2 \ ITT-b 



1 2 c L<>vap H >v 



(A.l) 



in terms of three invariant amplitudes Si, S 2 and S3, 
Vckm is a CKM factor. The quantities Si, S 2 and S3 
may be complex and involve two types of phases, CP- 
conserving strong phases and CP-violating weak phases. 
In general, the invariant amplitudes are a sum of several 
interfering amplitudes, S±j, S 2 j and Saj, respectively. 



From the decomposition Eq. (|A.4|) one finds the fol- 
lowing relations between the helicity amplitudes and 
the invariant amplitudes Si, S 2 , S3: 



1 



2miTO2 



(1 



Pn 



H± = S 1 ±f^\(l,mlm 2 2 ), 



rh 2 ) Si 



S2., „ 2 - 2- 
- — \{l,m 1 ,m 2/ 

S3 



(A.5) 



with A(l, m 2 , m|) = (1 — rh 2 ) 2 — 2m|(l + rh 2 ) + m\ and 
™i(2) = m 1(2) /m B . 

Note that the polarized decay amplitudes can be 
expressed in several different but equivalent bases. For 
example, the helicity amplitudes can be related to the 
spin amplitudes in the transversity basis (Ao , An , Aj_) 
defined in terms of the linear polarization of the vector 



mesons via: 



Ao = H , A 



H+ + P_ 



V2 



A 



P_ 



n/2 
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Table 5 Branching ratio [H, and decay amplitudes for B° -» if* p° [H, B° -> If* a; [H and B° -> if* 0, -»■ if* J/-0, 
B° -> iC* V(2S) [5|| . 



10 6 BR(£° -> K*°V) 

K\ 2 

axg{hj/h%) (rad) 

arg(feX/^) (rad) 
10 4 |S^| 
W 4 \S%\ 

w 4 \s%\ 

5\ - 5% (rad) 
S% - S% (rad) 
6% - tfV (rad) 



p° 




<t> 


J/i> 


i>{2S) 


3.4 


2.0 


9.8 


1330 


610 


0.70 


0.75 


0.480 


0.570 


0.46 


0.14 


0.12 


0.24 


0.219 


0.30 


1.17 


1.79 


2.40 


-2.86 


-2.8 


1.17 


1.82 


2.39 


3.01 


2.8 


1.17 


0.81 


2.66 


33.64 


28.86 


2.65 


1.67 


5.20 


42.49 


52.65 


2.31 


1.64 


5.28 


115.28 


153.00 


1.17 


1.79 


2.40 


-2.86 


-2.8 


-2.11 


-1.53 


-0.84 


0.90 


1.62 


1.17 


1.82 


2.39 


3.01 


2.8 



A , Am, A_l are related to Si, S 2 and S 3 of Eq. (jA.ll) 



via 



A = --J—((l-m 2 1 -m 2 2 )S 1 
2mi m9 V 



2fn\iri2 
S 

-— A(l,ra?,m|)), 



A,, =V2 5i, Aj. 



A(l,mf,m|) 



5, 



(A.6) 



The amplitude A\ (A = 0, || ,_L) are related to the in- 
variant amplitudes of the B — > V\ Vi decay by the for- 
mulas 

A Q = -—^((l-m 2 1 -m 2 2 )S 1 
2mim2 V 

St 



A {1 =V2S U il = _^ ^4M g 3 , (A . 7) 

If the B — > V\ V2 decay is invariant under CP transfor- 
mation, then Ao = Aq, A|| = A||, and Aj_ = — Aj_. 
The decay width is expresses as follows: 



r(B - Vl V 2 ) = -WA(l,m^i) |ycKM| 

x(|A | 2 + |A||| 2 + |A ± | 2 ) . (A.8) 

The matrix element for the — > K*° V decay, 
where V = p° ,u) ,<f> , J/i/)(lS) ,tp(2S) , . . . mesons, we 
can represent as 

V 777 p 



M = 



V2 



(A.9) 



Next, we define the normalized amplitudes: 
A A 



\/£v MM 2 
J>a| 2 = 1 (A,A' = 0, 



(A.10) 



By putting mi = my, 7712 = and using (|A.8|) . 

(jA.lOp we obtain the relation between the amplitudes 
h\ and A\ of the process under study B® — > K*° V for 
any vector meson V = p° , ui , <p , J/ipO-S) , tp(2S) , . . .: 



h v = ^L ]VtlVts 



4V2 '"'y irBR(B° d -> K*°V) 
x\^(l,m 2 v ,m 2 K «)AY, 



(ATI) 



where BR(. . .) is the branching ratio of B® —> K*° V 
decay and tb is the lifetime of a B meson. 

Solving Eqs. (|A~6) we find the scalars S\ , S2 and 6*3, 
and then extend the helicity amplitudes A\ off the mass 
shell of the meson V, i.e. for q 2 ^ rriy. We introduce the 
phases 8\ = axg{hj), SY = arg(Sf ), where i = 1 , 2 ,3. 
Then we have 



1^1 



141 



l#l 



V2 



A(l,m 2 / ,m| f ,) 
+ (l - my — m K 
+4V2rriK* my (1 



V A(l,m 2 / ,m 2 - t ) 
Sm^rn^A^ 2 



I All, 



- 2 ~2 

m v - m K , 
1/2 



sin((5 2 / - 



COs(SY — £q 



x\AY\\AJ\co S (6J-8X) 
V2 

x (1 - m 2 v - m 2 K ,)\Aj\ sin( ^r ~ S o) ' 
V2 



x((l-m 2 ,-m 2 f O|A[|cos( ( 5[-^) 
+2V2m y m K *|A^|) , 
^ = 5j (mod2ir),6Y = 5\ (mod2ir). (A.12) 
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